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CONSTITUTIVE  RELATIONS  FOR  AN  INELASTIC  GRANULAR  MEDIUM^ 


By  J.  B.  Weidler,  Jr.^  and  P.  R.  Paslay^ 

ABSTRACT 

The  results  of  an  experimental  and  analytical  program  to  develop  stress- 
strain  relations  for  a  granular  medium  are  presented.  Explicit  three-dimen¬ 
sional  relations  are  given  which  represent  the  inelastic  deformation  of  the 
medium.  A  simple  version  of  the  formulation  is  used  to  represent  experimental 
results  for  a  locally  quarried  quartz,  and  an  evaluation  of  the  relations  is 
presented.  The  simplified  relations  are  seen  to  be  a  compromise  between  math¬ 
ematical  tractability  on  the  one  hand  and  accuracy  of  prediction  on  the  other. 
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CONSTITUTIVE  RELATIONS  FOR  AN  INELASTIC  GRANULAR  MEDIUM^ 

By  J.  B.  Weidler,  Jr.^  and  P.  R.  Paslay^ 

INTRODUCTION 

Technological  advances  in  soil  mechanics  have  followed  a  different  devel¬ 
opment  pattern  from  strength  of  materials  and  elasticity  theory  for  structural 
and  machine  design  applications.  While  elasticity  theory  is  capable  of  predict¬ 
ing  most  structural  and  machine  design  behavior  within  design  accuracy,  there 
is  no  corresponding  theory  for  soil  mechanics. 

Mathematical  description  of  soil  behavior  is  extremely  difficult.  In  addi¬ 
tion  to  having  elastic  characteristics  under  certain  loadings,  its  behavior  can 
be  plastic,  viscoelastic  or  thixotropic  under  other  situations.  The  development 
of  a  mathematical  theory  which  encompasses  all  these  characteristics,  even  as  an 
engineering  approximation,  has  not  been  accomplished.  The  acceptable  mathemati¬ 
cal  complexity  of  any  material  description  is  limited  by  the  ability  to  solve 
practical  boundary  value  problems  with  the  description.  Although  numerical  meth 
ods  and  digital  computers  have  recently  extended  the  acceptable  mathematical  com 
plexity,  there  are  considerations  in  soil  mechanics,  such  as  material  stability, 
that  are  still  beyond  these  extended  capabilities. 

Most  efforts  to  deduce  useful  mathematical  results  for  practical  soil  me¬ 
chanics  problems  fall  into  two  categories.  The  first  of  these  is  elementary, 
three  dimensional,  theories  which  describe  a  limited  type  of  behavior.  The 
idealizations  of  Elasticity,  Plasticity  and  Coulomb  theory  represent  examples 


Work  presented  in  this  paper  was  sponsored  by  ARPA  Contract  SD-86. 
Brown  and  Root  Inc.,  Houston,  Texas,  formerly  Brown  University. 
Brown  University,  Providence,  Rhode  Island. 
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of  this  approach.  The  results  obtained  using  these  methods  in  soil  mechanics 
have  been  more  helpful  in  establishing  trends  rather  than  providing  quantita¬ 
tive  answers . 

The  second  category  is  equations  which  have  been  deduced  empirically  to 
predict  behavior  under  certain  circiomstances .  Most  design  formulas  fall  into 
this  category.  The  results  in  this  category  often  have  a  theoretical  basis 
but,  in  the  long  run,  rely  on  empiricism.  This  approach  in  soil  mechanics  has 
provided  the  technical  basis  for  practically  all  design  and  will  probably  con¬ 
tinue  to  play  this  role  for  many  years. 

The  work  presented  here  falls  into  the  first  category.  The  idealization 
under  consideration  is,  broadly  speaking,  a  granular  medium.  This  means  the 
behavior  is  time  independent  (the  stress-strain  equations  do  not  have  a  char¬ 
acteristic  time).  A  pore  pressure  is  not  included  in  this  formulation  although 
it  is  thought  to  be  straightforward  to  include  this  phenomenon  in  the  theory. 

Even  in  this  restricted  area  adequate  mathematical  representations  of  ma¬ 
terial  behavior  is  lacking.  The  most  common  representation  for  granular  mate¬ 
rials  is  the  Coulomb  theory.  The  basis  for  this  theory  is  that  whenever  on  any 
plane  the  ratio  of  resolved  shearing  stress  to  the  sum  of  the  normal  stress  and 
a  constant  called  the  "cohesive  strength"  reaches  a  critical  value ,  the  mate¬ 
rial  flows.  The  description  of  how  the  material  flows  is  the  part  of  this  the¬ 
ory  that  is  inadequately  resolved.  One  approach  which  is  used,  for  example,  in 
Terzaghi's  text^  is  to  assume  the  flow  is  incompressible.  It  is  well  known  that 
granular  media  can  expand  or  compact  during  flow.  Another  objection  to  the  in¬ 
compressibility  assumption  is  that  Drucker^  has  shown  that  normality  of  the 

^  Terzaghi,  K. ,  Theoretical  Soil  Mechanics ,  Wiley,  New  York,  1943. 

^  Drucker,  D.  C.,  "A  More  Fundamental  Approach  to  Plastic  Stress-Strain  Relaticns," 
Proc.  First  U.  S.  Nat.  Congress  of  App.  Mech.,  ASME,  1951,  pp.  487-491. 
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deformation  rate  vector  to  the  loading  surface  is  required  if  material  stabil¬ 
ity  is  to  be  insured.  The  condition  of  normality  necessitates  a  volume  change 
and  this  is  the  basis  for  a  second  approach  to  the  description  of  the  material 
flow.  A  normality  requirement  leads  to  a  plasticity  type  formulation  for 
soils 3.  Hansen**  has  conducted  experiments  to  try  to  evaluate  the  normality 
assumption  for  soils.  He  observed  volume  changes  which  were  between  the  zero 
volume  change  required  by  an  incompressibility  assumption  and  the  change  re¬ 
quired  for  normality.  The  significance  of  these  results  are  that  1)  there 
are  appreciable  volume  changes  during  flow  and  2)  unstable  material  behavior 

should  be  anticipated  for  certain  flow  conditions. 

5,6 

Earlier  work  by  the  authors  on  granular  media  was  directed  to  trying  to 
develop  a  plasticity  type  theory  with  a  modified  "stability"  requirement  so 
that  normality  was  not  a  consequence  of  the  requirement.  The  reason  such  an 
approach  is  attractive  is  that  the  stability  postulate  due  to  Drucker  has  a 
physical  basis  which  can  be  deduced  firom  energy  considerations.  In  the  authors’ 
earlier  work  it  was  postulated  that  energy  can  be  stored  and  released  in  the  ma¬ 
terial.  The  energy  is  taken  to  be  a  function  of  the  specific  volume.  The  re¬ 
sult  of  this  study  was  that  the  amount  of  energy  storage  and  release  would  be 
greater,  in  order  to  explain  Haythomthwaite ' s  results,  than  one  would  expect 
on  physical  grounds. 


Drucker,  D.  C.,  Henkel,  D.  J.  and  Gibson,  R.  E.,  "Soil  Mechanics  and  Work- 
Hardening  Theories  of  Plasticity,"  ASCE  Trans . ,  Vol.  122,  1957,  pp.  338-346. 

Hansen,  J.  B.,  Earth  Pressure  Calculation,  Danish  Technical  Press,  Copen¬ 
hagen,  Denmark,  1953. 

Weidler,  J.  B.  and  Paslay,  p.  R.,  "Analytical  Description  of  Behavior  of 
Granular  Media,"  Jour,  of  Engr.  Mech.  Div.,  Proc.  ASCE,  April  1969,  Vol.  95, 
EM2,  pp.  379-395. 

Paslay,  P.  R.  and  Weidler,  J.  B.,  "An  Analysis  of  the  Triaxial  Test  for  Co¬ 
hesionless  Soils,"  Jour,  of  Engr.  Mech.  Div.,  Proc.  ASCE,  June  1969,  Vol.  95, 
EM3. 


-4- 


The  work  presented  here  is  a  different  attempt  to  formulate  stress-strain 
equations  for  granular  materials  which  are  predictive  and  yet  tractable.  The 
notion  of  a  limit  surface  has  been  modified  to  make  the  presentation  of  this 
theory  straightforward.  A  relationship  is  introduced  between  the  mean  normal 
stress,  the  density  and  the  rate  of  change  of  volume.  The  result  is  a  specific 
three-dimensional  theory  that  predicts  several  phenomena  which  earlier  theories 
could  not.  Experimental  results  from  shear  tests  are  used  to  give  a  critical 
evaluation  of  the  theory  which  suggests  it  is  an  engineering  compromise  between 
accuracy  and  complication  which  should  be  adequate  for  many  practical  situations. 

THEORY;  GENERAL  CONSIDERATIONS 

This  development  depends  on  several  physical  assumptions  which  are  stated 
here  at  the  onset  so  the  mathematical  development  can  proceed  without  interrup¬ 
tion  . 

When  an  element  of  a  granular  medium,  large  enough  to  be  treated  as  a  con¬ 
tinuum,  is  loaded  monotonically ,  it  responds  elastically  and  it  is  assumed  that 
one  of  two  phenomena  occur  when  the  load  becomes  sufficiently  high.  The  first 
is  that  the  individual  grains  crush.  This  will  occur,  for  example,  if  the  load¬ 
ing  is  hydrostatic  compression.  The  second  is  that  the  element  distorts  due  to 
shearing  motion  as  in  a  shear  test.  Undoubtedly  some  particles  fracture  during 
shearing,  but  this  is  neglected.  The  only  response  which  this  theory  considers 
is  the  shearing  response. 

The  magnitude  oF  the  shear  stress  required  to  initiate  flow  is  assumed  to 
be  a  function  of  the  mechanical  state  of  the  material  and  mean  normal  stress. 

This  mechanical  state  is  assumed  to  be  measured  by  the  specific  volume  of  the 
material . 
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In  conventional  plasticity  theory  the  yield  surfaces  are  all  cylinders  in 
principal  stress  space  and  are  identified  by  their  cross  section  on  the  ir  plane 
(the  Mises  condition  is  a  circle  while  the  Tresca  condition  is  a  hexagon ) .  The 
Coulomb  condition  is  a  six  sided  pyramid  in  principal  stress  space.  On  the  ir 
plane  a  six-sided  figure  with  threefold  symmetry  results.  In  the  theory  devel¬ 
oped  here,  a  surface  of  revolution  will  represent  the  limit  surface  at  any  time. 
The  intersection  of  a  ir  plane  with  the  surface  of  revolution  will  be  a  circle 
whose  radius  depends  on  the  mechanical  state.  This  choice  of  limit  surface 
makes  possible  an  approximation  to  Coulomb  theory  and  is  mathematically  easier 
to  work  with  than  the  Coulomb  theory  for  many  problems. 

Experiments ,  such  as  shear  box  tests ,  indicate  that  when  flow  is  induced , 
there  is  a  volume  change.  As  shearing  proceeds,  the  volume  reaches  an  equilib¬ 
rium  value  which  is  related  to  the  normal  load  applied  during  shearing.  If  the 
specific  volume  is  above  (below)  its  equilibrium  value  when  shearing  commences, 
then  the  specific  volume  decreases  (increases)  during  the  shearing  flow  and  ap¬ 
proaches  the  equilibrium  value. 

In  order  to  incorporate  the  above  ideas  into  a  theory,  it  is  convenient  to 
work  in  a  strain  space  rather  than  a  stress  space  owing  to  the  mechanical  state 
being  described  in  terms  of  specific  volume.  The  reduced  elastic  components  of 
strain  are  deduced  from  the  difference  of  the  total  reduced  strain  e  and  the 
reduced  inelastic  strain  a  .  The  stresses  and  elastic  strains  are  to  be  rela¬ 
ted  through  a  linear  Hooke's  law  so  that  the  Mises  condition  may  be  written  as 

(e  -  a)*(e  -  a)  $  p2(a,o) . (1) 

where  p(o, a)  abbreviated  as  p  is  the  elastic  shear  strain  required  to  initi¬ 
ate  inelastic  flow,  a  is  the  mean  normal  stress  (tensile  stresses  measured  as 


-6- 


positive)  and  a  is  the  inelastic  volumetric  strain.  The  function  p(a,a) 
will  be  assumed  to  have  the  following  properties 


p(0,a)  =  0 . (2) 

P(o,0)  =  0 . (3) 


The  first  property  simply  implies  that  when  there  is  no  mean  normal  compres¬ 
sion,  there  is  no  resistance  to  flow.  The  second  property  adjusts  the  refer¬ 
ence  state  from  which  the  plastic  volume  change  is  measured.  Fig.  1  depicts 
the  loading  surface  graphically. 

A  second  function  is  introduced  to  determine  the  equilibrium  volume  for 
a  given  mean  normal  stress.  The  function  is  denoted  by  a°(a)  or  abbrevia¬ 
ted  as  a°  . 

The  requirements  for  inelastic  flow  to  occur  are  that  the  Mises  condition 
equality 

(e  -  a)*(e  -  a)  =  p2(o,a) .  (4) 

be  satisfied  and  that  a  loading  condition,  deduced  by  differentiation  of  the 
Mises  condition,  inequality, 

(e  -  a)*e  -  j  a  >  0 . (5) 

be  satisfied.  If  Eqs.  4  and  5  are  not  simultaneously  satisfied,  then  only 
elastic  deformation  occurs. 

When  inelastic  flow  does  occur,  a  flow  law  is  required  to  describe  the 
deformation.  To  describe  the  flow  law  Introduced  here,  it  is  convenient  to 
refer  to  a  principal  strain  space  where  the  Mises  condition  describes  a  cyl¬ 
inder.  If  the  loading  conditions,  Eqs.  4  and  5,  are  satisfied,  then  the  di¬ 
rection  of  the  deformation  rate  vector  in  the  strain  space  is  such  that  its 
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projection  on  the  tr  plane  going  through  point  A  is  normal  to  the  cylinder. 
The  component  of  the  deformation  rate  vector  perpendicular  to  the  ir  plane  is 
determined  by  the  mean  normal  stress  and  the  separation  between  the  current  in¬ 
elastic  volume  strain,  a  ,  and  the  equilibrium  volume  strain  a°(a)  .  Analy¬ 
tically  this  may  be  expressed  as  follows 

a  =  X(e  -  a)  if  yield  and  loading  4  “  ®  . 

othewise 

a  =  Xh(a,a®  -  a)  conditions  are  satisfied,  a  =  0  . (7) 

In  Eqs.  6  and  7,  X  is  an  arbitrary  parameter  and  h(a,a°  -  a)  abbreviated 
as  h  is  a  function  which  has  the  following  property 

h(a,0)  =  0 . (8) 


since  the  inelastic  volume  change  and  the  equilibrium  volume  change  coincide 
in  this  circumstance.  X  may  be  eliminated  from  Eqs.  6  and  7  by  differenti¬ 
ating  the  Mises  condition  equality  (which  must  remain  valid  during  inelastic 
flow)  and  substituting  Eqs.  6  and  7.  This  procedure  leads  to 


X 


(e-a)«e  -  j  (p^)*^  o 
p2  +  i  (p2),^  h 


(9) 


so  that,  by  noting  Eq.  5,  X  is  positive  when  p2  +  ^  ^  positive. 

When  X  is  positive,  it  simply  means  that  the  projection  of  the  deformation 
rate  vector  on  the  ir  plane  at  point  A  is  directed  radially  outward  from  the 
Mises  condition  surface.  Substitution  of  Eq.  9  into  6  and  7  gives 
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(e-a)*e  -  ^  ° 

- - - —  (e  -  a) 

(e-a)*e  -  i  (P^)»„  o 
— : — :: — £ - h 

p2  +  I  (p2),^  h 


if  yield  and 
loading  con¬ 
ditions  are 
satisfied 


a  =  0 


(10) 


otherwise 


a  =  0 


(11) 


Finally,  the  stress  and  elastic  strain  states  must  be  related  through 
Hooke's  law.  The  reduced  elastic  strain  is  (e  -  a)  and  the  elastic  volume 
change  per  unit  volume  is  (e  -  a)  so  that 


s  =  2G(e  -  a) . (12) 

o  =  K(e  -  o) . (13) 


where  K  and  G  are  functions  of  a  and  o  . 

Eqs.  4,  5,  10,  11,  12  and  13  complete  the  constitutive  formulation  for 
the  material.  For  the  solution  of  boundary  value  problems,  the  equilibrium 
equations  and  strain  compatibility  equations  must  be  included  to  form  an  ap¬ 
propriate  set  of  field  equations. 


THEORY;  SIMPLEST  CASE 

In  order  to  make  a  practical  evaluation  of  the  formulation  presented  in 
the  preceding  section ,  a  specific  choice  has  been  made  for  each  of  the  func¬ 
tions  p(a,o),  o®(o)  and  h(a,a®  -  o)  .  The  simplest  case  which  yields 

qualitatively  correct  predictions  for  elementary  deformations  is  as  follows 


G(a,a)p(a,a)  =  +aoG  p  . (14) 

o  11 

a°(0)  =  a®  +  a®a  .......  . . (15) 

o  1 
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h(o,a®  -  a)  =  (a®  -  a)h^ . (16) 

where  >  “o’  “l  constants  dependent  on  the  particular 

granular  medium  being  represented.  The  Eqs.  14,  15,  16  satisfy  the  Condi¬ 


tions  2,  3,  8. 

With  the  simplified  definition  loading  occurs  whenever 

G^(o,a)(e  -  a)  (e  -  a)  =  G  . (17) 

-  ^  ^  ^  O  XX 

and 

oG,  * 

G2(a,a)(e  -  a)*e  -  G  - — -^)  ^  >  0 . (18) 


otherwise  there  is  no  loading.  The  stress -deformation  relations  then  become 


013  ,  • 

.  (G2(a,a)(e-a).e/G3,  2„2o2)  +  (—2--  1)^ 

a  __  ( 

1  +  h  (a®+a?a-a)— (1  -  -p --) 

1  o  1  a  G 


e  -  a) . 


oG,  • 

^  (G2(a,a)(e-a)*e/G^2p^^2„2jj2)  _  (i - - — )_ 


loading 


-  IXO, 

1  +  hda°+a?a-a)— (1 - 

1  o  1  a  G 


h  (a®  +  a?a  -  a) 
1  O  X 


a  =  0  . (21) 

)  no  loading 

i  =  0  . (22) 

In  the  following  sections  these  stress-deformation  relations  are  evalua¬ 
ted  against  experimental  results .  The  form  chosen  for  G(a ,a )  is  one  which 
has  been  determined  experimentally,  reported  in  the  literature^  and  can  be  ap¬ 
proximated  by 

^  Hardin,  Bobby  0.,  "Dynamic  Versus  Static  Shear  Modulus  for  Dry  Sand,"  ASTM 
Materials  Research  and  Standards,  Vol.  5,  No.  5,  pp.  232-235. 
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G(a,a)  =  10,850  (.0833  -  a)  ,  o  <  0 . (23) 


where  o  and  G(o,a)  are  in  psi. 


COMPARISON  OF  SIMPLEST  THEORY  WITH  EXPERIMENTAL  RESULTS 


The  comparison  with  experiment  to  be  presented  here  is  based  on  77  shear 
tests  which  were  run  on  locally  quarried  quartz  sand.  The  shear  test  apparatus 
is  shown  graphically  in  Fig.  2.  The  tests  used  in  this  evaluation  are  part  of 
a  large  set  of  tests  run  as  part  of  another  investigation  and  have  been  re¬ 
ported  in  detail® .  The  test  consists  of  shearing  the  top  half  of  a  shear  box 
at  a  slow  enough  rate  to  make  inertia  effects  negligible.  The  shearing  force, 

T  ,  required  to  accomplish  this  motion  is  measured.  The  normal  load  N  is  ap¬ 
plied  through  weights  so  that  it  remains  constant  during  the  test.  The  verti-^ 
cal  and  horizontal  components  of  the  shear  box  motion  are  measured. 

There  are  many  drawbacks  to  using  such  tests  for  the  evaluation  of  consti¬ 
tutive  relations.  The  most  obvious  is  that  the  state  of  deformation  in  the 
sample  is  not  uniform.  It  will  be  observed  in  the  following  that  a  number  of 
assumptions  is  required  in  order  to  reach  a  critical  evaluation  of  the  theory 
from  these  tests. 

A  typical  set  of  data  for  a  compacted  sand  is  sketched  in  Fig.  3a.  If  the 
data  beyond  about  two  grain  diameters  is  observed,  it  is  seen  that  the  average 
specific  volume  decreases  while  the  shear  load  also  decreases .  As  the  shearing 
progresses,  both  T  and  6^  appear  to  as5rmptotically  approach  constant  values. 
The  authors  have  ignored  6  displacements  of  less  than  two  grain  diameters  be- 
cause  extraneous  settling  and  orientation  effects  are  occurring  which  would 

8 

Weidler,  J.  B.,  "The  Flow  Behavior  of  Sand  at  Failure,"  Brown  University,  Di¬ 
vision  of  Engineering  Report  ARPA  E59,  August  1968. 
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possibly  not  be  present  in  a  natural,  undisturbed  foundation.  When  a  loose 
sand  is  tested,  the  results  follow  the  pattern  shown  in  Fig.  3b.  For  each 
test  the  following  data  were  recorded: 

initial  specific  volume,  v  (in^/#) 
nominal  normal  stress,  N  t  A 

nominal  initial  shearing  stress,  T  v  A  or  T.  ,.  t  A  (^/in^) 

max  2  gr*  dia 

2  gr.  dia 

The  experiments  reported  here  were  run  for  5  values  of  N  t  A  which  were  46.1, 
37.3,  28.5,  19.7  and  10.9  psi  and  for  assorted  values  of  v  from  17.5  to 
20.5  in^/lb.  It  should  be  noted  that  out  of  the  77  tests,  only  one  had  the 
character  of  Fig.  3b  owing  to  the  difficulty  of  obtaining  sufficiently  loose 
initial  states . 

Fig.  4  shows  the  result  of  plotting  the  nominal  initial  shearing  stress 
as  a  function  of  initial  specific  volume  for  the  different  values  of  nominal 
normal  stress.  An  approximation  to  the  value  of  G(a)pj^^  can  be  achieved  from 
this  curve  if  it  is  assumed  that  the  nominal  initial  shearing  stress  defines 
the  initiation  of  inelastic  flow  for  a  mean  normal  stress  equal  to  the  nominal 
normal  stress  and  a  volumetric  inelastic  strain  given  by 

a  =  ^^2^  units  of  in^/#) . (24) 

Eq.  24  is  deduced  from  the  figure  by  1)  extrapolating  the  curves  for  differ¬ 
ent  values  of  N  t  A  in  Fig.  4  to  a  nominal  initial  shearing  stress  of  zero, 

2)  finding  that  the  average  v  determined  this  way  for  the  different  values 
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of  N  V  A  is  24  in  ^/Ib  and  3 )  measuring  the  volumetric  strain  from  the  ref¬ 
erence  value  of  24  in^/lb.  The  above  scheme  then  leads  to  the  following  equa¬ 
tion: 

nominal  initial  shearing  stress  =  ^  A)(-o).  ....  (25) 

The  straight  lines  shown  in  Fig.  4  are  deduced  by  setting 

G  p,  =  4.00  . (26) 

O  11 

The  evaluation  of  aP  and  a°  ,  in  Eq.  15  can  be  accomplished  by  plotting 

dfiy  °  . 

3-t—  as  a  function  of  a  for  each  different  nominal  normal  stress.  In  this 
case  five  curves  are  obtained .  By  extrapolating  these  curves  to  the  point  where 

d6y 

is  zero,  an  approximation  to  the  equilibrium  volume  for  each  nominal  normal 
stress  may  be  obtained.  The  necessary  data  for  accomplishing  this  procedure  are 
plotted  in  Fig.  5.  A  plot  of  the  equilibrium  volumes  as  a  function  of  nominal 
normal  stress  which  results  from  this  procedure  is  given  in  Fig.  6.  If  N  A 
is  interpreted  as  -o  ,  the  straight  line  approximation  shown  in  Fig.  6  is  rep¬ 
resented  by 

a®  =  -  .124  . (27) 

o 

a®  =  +.0009  1/psi . (28) 

dSy 

To  make  use  of  measured  values  of  xr"  for  the  determination  of  G(o,a)h,, 

d6H  1 

a  brief  derivation  is  required.  For  the  first  increment  of  flow  from  Some 
specified  state ,  the  stress-deformation  Relations  19  and  20  predict  the  rates 
of  change  of  strain  and  volume.  If  the  ratio  of  the  xy  component  of  a  from 
Eq.  19  to  o  from  Eq.  20  is  formed,  the  following  equation  results 


-13- 


a  (e  -a  ) 

xy  _  _ xy  xy 

a  h  (o°+a?o-a) 

1  o  1 


This  equation  then  determines  G(0,a)h^  for  each  shear  test  in  terms  of  the 
constants  »  ®o  *  °i  have  already  been  determined  and  the  meas¬ 

ured  values  o  ,  o  and  d6,,/d6i,  .  As  before  a  is  approximated  by  -N  t  A 
and  the  procedure  is  straightforward.  The  value  given  for  G(o,a)  from  Eq.  23 
then  makes  possible  a  determination  of  h^^  for  each  experiment.  As  h^  was 
supposed  constant  in  the  theory,  this  calculation  procedure  furnishes  a  crucial 
test  for  the  theory.  Fig.  7  shows  the  results  of  the  calculations  for  h^^  . 

Whether  or  not  Fig.  7  is  regarded  as  a  partial  verification  of  the  theory 
or  not  is  a  matter  of  viewpoint.  In  classical  plasticity  theory  the  predicted 
and  measured  values  of  strain  increment  ratios^  vary  by  as  much  as  2:1  and  yet 


Naghdi,  P.  M.,  Rowley,  J.  C.,  and  Beadle,  C.  W.,  "Experiments  Concerning  the 
Yield  Surface  and  the  Assumption  of  Linearity  in  the  Plastic  Stress-Strain 
Relations,"  Journal  of  Applied  Mechanics,  Vol.  22,  No.  3,  1955,  pp.  416-420. 
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the  theory  is  widely  accepted.  It  appears  that  deviations  of  the  same  order 
are  present  here  for  a  much  more  complicated  material  behavior  so  that  the 
authors  view  Fig.  7  as  encouraging. 
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APPENDIX  -  NOTATION 


A  =  cross-sectional  area  of  shear  test  apparatus; 

see  Fig.  2 

a  =  reduced  inelastic  strain  tensor 

=  quadratic  invariant  of  b,  b*b=J^b..b.. 

i  j 

where  b^^  are  components  of  b  referred  to 
Cartesian  coordinates 


e  =  reduced  total  strain  tensor 

G  =  G(a,0)  =  shear  modulus  which  depends  on  inelastic  volume 

change  and  mean  normal  stress;  see  Eq.  12 

G^,  a®,  aj,  hj^  =  a  set  of  constants  which  describe  the  simplified 

theory.  See  Eqs.  1,  2,  3  and  10 


h  =  h(o,  a®  -  a) 
K  =  K(o,a) 

N 

s 

T 

V 

a 

a®  =  o®(o) 

e 


=  a  function  first  appearing  in  Eq.  7 

=  bulk  modulus  which  depends  on  inelastic  volume 
change  and  mean  normal  stress;  see  Eq.  13 

=  normal  force  on  shear  test  apparatus;  see  Fig.  2 

=  reduced  stress  tensor,  s..  =  a.  .  -  (iT  a,  ,)6.. 
v  -  •  *1-  ^  1]  1]  3  f  kk  in 

where  a  is  the  stress  k 

=  tangential  force  on  shear  test  apparatus ;  see 
Fig.  2 

=  specific  volume 

=  inelastic  volume  strain  measured  from  loose  state 

=  a  function  first  appearing  in  Eq.  7 

=  displacement  downward  of  shear  test  apparatus  cover 
see  Fig .  2 

=  horizontal  displacement  of  shear  test  apparatus 
cover;  see  Fig.  2 

=  total  volume  change  measured  from  loose  state 

=  parameter  first  appearing  in  Eqs.  6  and  7 
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p  =  p(a,a)  =  function  of  mean  normal  stress  and  inelastic 

volume  change  first  appearing  in  Eq.  1 

o  =  mean  normal  stress 


CAPTIONS  FOR  FIGURES 


Fig.  1 

Fig.  2 
Fig.  3 

Fig.  4 

Fig.  5 

Fig.  6 

Fig.  7 


-  Sketch  showing  relation  of  strain  point  to  limit  surface  in 
strain  space. 

-  Shear  box  testing  apparatus . 

-  Sketches  of  Typical  Results  from  Shear  Box  Tests,  a)  for  a 
Compacted  Sand  and  b)  for  a  Loose  Sand. 

-  Nominal  initial  shearing  stress  versus  initial  specific  volume 
for  indicated  values  of  the  nominal  normal  stress. 

-  Rate  of  change  of  the  vertical  displacement  with  respect  to  the 
horizontal  displacement  as  a  function  of  the  initial  inelastic 
volume  strain  measured  from  loose  state  for  indicated  values  of 
the  nominal  normal  stress. 

-  Equilibrium  inelastic  volume  strain  as  a  function  of  the  nominal 
normal  stress . 

-  Variation  of  h,  with  inelastic  volume  strain  and  nominal  normal 

1 

stress.  For  a  perfect  match  between  theory  and  experiment  h^^ 
is  constant. 
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